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ABSTRACT 

We have studied the globular cluster systems of 7 giant, edge-on spiral galaxies using Hubble 
Space Telescope imaging in V and /. The galaxy sample covers the Hubble types Sa to Sc, al- 
lowing us to study the variation of the properties of globular cluster systems along the Hubble 
sequence. The photometry reaches ^1.5 mag beyond the turn-over magnitude of the globu- 
lar cluster luminosity function for each galaxy. Specific frequencies of globular clusters {Sn 
values) were evaluated by comparing the numbers of globular clusters found in our WFPC2 
pointings with those in our Milky Way which would be detected in the same spatial region if 
placed at the distance of the target galaxies. Results from this method were found to be con- 
sistent with the more commonly used method of constructing radial distribution functions of 
globular clusters. The Sn values of spirals with B/T ^ 0.3 (i.e., spirals with a Hubble type 
later than about Sb) are consistent with a value of Sm = 0.55±0.25. We suggest that this pop- 
ulation of globular clusters represents a 'universal', old halo population that is present around 
each galaxy. Most galaxies in our sample have Sn values that are consistent with a scenario 
in which globular cluster systems are made up of ( i) the aforementioned halo population plus 
( ii) a population that is associated with bulges, which grows ^ linearly with the mass of the 
bulge. Such scenarios include the 'merger scenario' for the formation of elliptical galaxies as 
well as the 'multi-phase collapse' scenario, but it seems inconsistent with the 'secular evolu- 
tion' scenario of Pfenniger & Norman (1990), in which bulges are formed from disc stars by 
means of the redistribution of angular momentum through bar instabilities and/or minor per- 
turbations. However, there is one bulge-dominated spiral galaxy in our sample (NGC 7814) 
with a low Sn value that is consistent with those of the latest-type spirals. This means that 
the 'secular evolution' scenario can still be viable for some bulge-dominated spirals. Thus, 
our results suggest that the formation histories of galaxy bulges of early-type spirals can be 
significantly different from one galaxy to another. 
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1 INTRODUCTION 

Recent observations with the Hubble Space Telescope (HST) and 
large-field ground-based CCD cameras have caused rapid advances 
in our knowledge of the formation and evolution of globular clus- 
ter (GC) systems of galaxies. However, ellipticals and SO galaxies 
have received by far the largest amount of attention, mainly because 
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early-type galaxies usually have much richer GC systems and also 
suffer less from internal extinction problems than do spirals. Conse- 
quently, our knowledge of the GC systems in spiral galaxies is still 
limited to a handful of galaxies (e.g., Harris 1991; Ashman & Zepf 
1998; Kissler-Patig et al. 1999). However, it should be recognized 
that this scarcity of data for spiral galaxy GC systems constitutes 
an important limitation to the use of GC systems as probes of the 
formation and evolution of both early-type and late-type galaxies. 

For instance, the ability to test predictions of scenarios for the 
formation of galaxies and their GCs depends on our knowledge 
of the 'typical' properties of GC systems of spirals. A good ex- 
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ample is the 'merger model' (Schweizer 1987; Ashman & Zepf 
1992), in which elliptical galaxies are formed by mergers of spi- 
ral galaxies. In this picture, the GC systems of elliptical galaxies 
are composite systems. One population of GCs is associated with 
the progenitor spirals (i.e., metal-poor GCs with a halo-like spa- 
tial distribution), while a second population of clusters forms in the 
merger event (i.e., metal-rich GCs with a bulge-like spatial distri- 
bution). Ashman and Zepf (1992; 1998) described several testable 
predictions arising from this scenario, most of which are however 
based on comparisons with properties of 'typical' spiral galaxy GC 
systems. A particular problem in this respect is that the character- 
istic specific frequency of GCs (i.e., the total number of GCs per 
unit galaxy luminosity) around spirals is poorly known, especially 
for late-type spirals (e.g., the compilation of Ashman & Zepf 1998 
includes specific frequency measurements for only 5 spirals with 
Hubble type Sb or later). The absence of firm constraints on the 
number of GCs contributed by the progenitor spirals leads to un- 
certainties in many of the key predictions for elliptical galaxy GC 
systems that follow from the merger model. Furthermore, metal- 
rich (—1.0 ^ [Fe/H] ^ —0.2) GCs are known to be associated 
with bulges of spiral galaxies like the Milky Way and M3 1 rather 
than with their discs or halos (Minniti 1995; Barmby, Holland & 
Huchra 2001; Forbes, Brodie & Larsen 2001), so that mergers of 
spirals containing significant bulges are likely to contribute metal- 
rich GCs to the merger remnant which are not formed during the 
merger. It is therefore very important to establish — using a sig- 
nificant sample of spiral galaxies — whether or not the properties 
of the Milky Way (MW) GC system are typical of spiral galaxy 
GC systems as a class, and to study the relationship between bulge 
luminosity and the number of metal-rich GCs in spiral galaxies. 

Another important and timely area where knowledge of the 
properties of GC systems of spirals can yield significant progress is 
that of the formation and evolution of bulges of spiral galaxies. On 
the one hand, structural and dynamical properties of spiral bulges 
have long shown strong similarities with those of (low-luminosity) 
elliptical galaxies. Bulges follow the Fundamental Plane of ellip- 
tical galaxies (e.g.. Bender, Burstein & Faber 1992), while their 
internal dynamics are consistent with oblate, isotropic models, just 
like low-luminosity ellipticals (e.g., Davies et al. 1983). Extinction- 
corrected optical and near-IR colours of bulges have been shown to 
be very similar to those of elliptical galaxies in the Coma cluster 
(Peletier et al. 1999), and spectroscopic metallicities and [Mg/Fe] 
ratios in bulges are similar to those of ellipticals at a given bulge 
luminosity (Goudfrooij, Gorgas & Jablonka 1999). As ellipticals 
in rich clusters most likely formed at redshift z > 3 (e.g., Stan- 
ford, Eisenhardt & Dickinson 1997), these similarities argue for 
an early formation of bulges, in line with the original 'monolithic 
collapse' model by Eggen, Lynden-Bell & Sandage (1962). An al- 
ternative (and quite different) model is that bulges form from disc 
material through redistribution of angular momentum (Pfenniger 
&. Norman 1990). In this scenario, large amounts of gas are driven 
into the central region of the galaxy by a stellar bar and trigger in- 
tense star formation. If enough mass is accreted, the bar itself will 
dissolve and the resulting galaxy will reveal a bigger bulge than 
before bar formation; galaxies would thus evolve from late to ear- 
lier types along the Hubble sequence (see also Pfenniger, Combes 
& Martinet 1994). Note that this scenario does not involve GCs, 
in that only disc stars would contribute to the secular building of 
bulges. Hence, there would be no reason for a relation between the 
bulge-to-total luminosity ratio and the number of GCs per galaxy 
in this scenario. This provides a testable prediction, and is one of 
the key tests performed in the present study. 



This paper is built up as follows. Section|2|describes the sam- 
ple selection and the HST observations, while the data analysis is 
described in section |3] The various results are discussed in sec- 
tion|4] Finally, section |5| summarizes the main conclusions of this 
study. 



2 APPROACH 

2.1 Galaxy sample selection 

The sample consists of spiral galaxies in an edge-on configuration. 
We selected such galaxies in order to enable detection of GCs on 
both sides of the spiral discs and to minimize (spatially) the impact 
of dust absorption. In addition, the edge-on configuration allows 
one to assess radial number density distributions of the GC systems 
around these galaxies. Northern galaxies were selected from the 
Uppsala General Catalogue (UGC; Nilson 1973), while southern 
galaxies were selected from the Surface Photometry Catalogue of 
ESO-Uppsala Galaxies (ESO-LV; Lauberts & Valentijn 1989). Our 
two first-cut selection criteria were: 

• the inclination of the galaxy i > 80° , and 

• The Galactic latitude fe > 30°. 

Inclinations were determined following Guthrie (1992), assuming 
an intrinsic flattening go = {b/a)o = 0.11. Inclinations were then 
derived using Hubble's (1926) formula 

cos^i = (q^ - go)/(l - ?o). 

where q — b/a is the observed (catalogued) axis ratio. 

In order to get the highest observing efficiency out of the HST 
observations (see below), we further considered two main (and mu- 
tually counteracting) factors. On the one hand, a galaxy should be 
near enough to reach Mv = —6.0, which is ~ 1.5 magnitudes be- 
yond the peak in the Galactic GC luminosity function (hereafter 
GCLF; Harris 1996), in a few HST orbits using WFPC2. On the 
other hand, a galaxy should be distant enough for a significant frac- 
tion of its GC system to fit in the field of view of one or two WFPC2 
exposures. 

From the remaining list of candidate spirals, we selected 2-3 
galaxies of each main Hubble type (Sa, Sb, and Sc) to study the 
variation of the properties of GC systems as a function of Hubble 
type (i.e., of bulge-to-total luminosity ratio). Global properties of 
the final sample galaxies are given in TableQ 

2.2 Observations 

Observations were made with the Wide Field and Planetary Cam- 
era 2 (WFPC2) aboard HST as part of General Observer program 
6685. The data consist of multiple images through the F555W and 
F814W filters. Our program was supplemented by archival images 
of a few sample galaxies (NGC 4565 and NGC 4594), taken from 
other HST programs. The images of NGC 4565 were not reana- 
lyzed; all relevant data on GC candidates in that galaxy was taken 
from Kissler-Patig et al. (1999). The HST observations are listed in 
Table|2|together with the exposure times for each galaxy. For a few 
galaxies, a subset of the images were spatially offset by 0'.'5 from 
the others (corresponding to an approximately integer pixel shift in 
both PC and WF CCDs of the WFPC2) to enable a good correc- 
tion for hot pixels. The locations of the WFPC2 fields are shown 
in Fig.Q superposed onto grey-scale images from the Digital Sky 
Survey. 
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Table 1. Global properties of the sample galaxies. 
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Notes. From Lyon Exti'agalactic Database (LEDA;"http : / / leda . univ-lyonl . f r i. Corrected for both internal and foreground reddening. 
Galactic foreground extinction only. References are given in galaxy order. Derived from the above values. ^ Corrected value (see Section lT2l 2). ^ AH 

bulge radii are expressed in terms of the equivalent radius of an ellipse (see text). Assumed half-light radius for the bulge (see Section|^2|2). Equal to 
either r^^v or 1.678 Toxp, depending on whether the radial light profile is fit better with r^/^ or with exponential profile. REFERENCES: (1) NASA/IR\C 
Extragalactic Database (NED); (2) de Vaucouleurs et al. (1991; hereafter RC3); (3) Schlegel et al. (1998); (4) Willick et al. (1997; using Tully-Fisher relation 
(TP)); (5) Tully & Pierce (2000; using TP); (6) Perrarese et al. (2000; using Cepheids); (7) Jacoby, Ciardullo & Harris (1996; using a variety of methods); (8) 
Pord et al. (1996; using Planetary Nebula Luminosity Function); (9) Tonry et al. (2001; using Surface Brightness Fluctuations); (10) this paper; (11) Gavazzi 
et al. (2000); (12) Wu et al. (2002); (13) Baggett, Baggett & Anderson (1998); (14) Bothun, Harris & Hesser (1992). 




Figure 1. Images of the sample galaxies. The greyscale images were extracted from the Digital Sky Survey, and the superposed lines show the locations of the 
WFPC2 field(s) of view (where each CCD chip is outlined individually). The galaxy identifications and field sizes are indicated in a comer of each image. 
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2.3 Data reduction and photometry 

We used the IRAF^ package for data reduction. After standard 
pipeline processing and alignment of the images, we combined 
the images using the IRAF/STSDAS task CRREJ, thereby effec- 
tively removing both cosmic rays and hot pixels. We also trimmed 
each image to exclude the obscured regions near the pyramid edges 
of WFPC2, yielding 751 x 751 usable pixels per CCD. Prior to 
performing source photometry, the strongly varying galaxy back- 
ground was fitted and subtracted. The main reason for this is to 
minimize errors in the photometry due to any particular choice of 
object aperture and sky annulus. For the CCD chips covering a sig- 
nificant part of the galaxies' bulges or dusty discs, we applied a 
median filter with a 30 x 30 pixel kernel to approximate the galaxy 
background. For the other CCDs, the smooth background gradient 
was fitted by a bi-cubic spline fit. 

The photometry of star cluster candidates was carried out us- 
ing the DAOPHOT-II IStetson \ WJ \ package within IRAF. The ob- 
jects were selected by applying the DAOFIND task to an image pre- 
pared by dividing the data image in question by the square root of 
the appropriate model image (created as described in the previous 
paragraph). This procedure ensures uniform shot noise characteris- 
tics over the whole image. We adopted fairly tight shape constraints 
[—0.6 < roundness < 0.6; 0.2 < sharpness < 0.9] in order to ex- 
clude extended background galaxies and faint objects distorted by 
noise or any residual bad pixels. The detection threshold was set 
at 4 sigma above the residual background. Although DAOFIND re- 
turned with apparent point-like detections located within - or on the 
edge of ~ dust features in the galaxy discs, we decided to exclude 
those from further analysis due to the difficulty in judging the loca- 
tion of those sources relative to the dust features along the line of 
sight. 

PSF photometry was performed using the ALL STAR routine. 
For each chip, a PSF was constructed using, when possible, 10-15 
bright, isolated GCs in the field. Aperture corrections were calcu- 
lated using growth curve analysis on a subset of the PSF objects on 
each chip. One advantage of PSF photometry over the commonly- 
used two-pixel aperture photometry is that it avoids the systematic 
photometric offset due to underestimating aperture corrections for 
bright, resolved GCs. 

One of the main sources of uncertainty in WFPC2 photome- 
try (in general) is due to charge-transfer efficiency (CTE) problems 
of the CCDs for which correction recipes are available (Whitmore, 
Heyer & Casertano 1999; Dolphin 2000). We applied the Whitmore 
et al. formulae to our photometry. However, in the present case of 
GC photometry on a relatively high background (due to the diffuse 
galaxy light), the CTE correction typically affected the magnitudes 
by less than 0.02 mag and the colours by less than 0.01 mag, negli- 
gible for the purposes of this paper. 

The measured WFPC2 STMAG magnitudes F547M, F555W, 
and F814W were converted into Johnson-Cousins V and / mag- 
nitudes using colour terms as derived by Holtzman et al. 1 1995 i 
and Goudfrooij et al. <2001t . Specifically, the following conver- 
sions were used: 

V - F547M = (0.009 ± 0.002) 

+ (0.005 ± 0.002) (F547M - F814W) 
- (0.008 ± 0.002) (F547M - F814W)^ 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. U.S.A. 



Table 2. Summary of the HSTAVFPC2 images of the sample galaxies. 
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Notes, a: Number of exposures; h: 2 WFPC2 pointings (exp. times listed 
are per pointing); c: Integer-pixel dithered data. 



V - F555W = (0.643 ± 0.001) 

+ (0.004 ± 0.002) (F555W - F814W) 
- (0.015 ± 0.002) (F555W - F814W)^ 

/-F814W = (-1.266 ±0.001) 

+ (0.018 ± 0.003) (F555W - F814W) 
+ (0.016 ± 0.002) (F555W - F814W)^ 

Finally, the magnitudes of the point-like sources in the sam- 
ple galaxies were corrected for Galactic foreground extinction (cf. 
Table 0. The Av values were converted to Aj according to the 
Galactic extinction law of Rieke & Lebofsky <1985> . 

2.4 Cluster candidate selection 

To select candidate GCs, we adopted a fairly generous selection cri- 
terion for the colour cut: 0.3 < V^ — / < 2.0. The lower limit was 
set to 0.2 mag bluer than the bluest GC in the Milky Way, while 
the upper limit was set to prevent the exclusion of somewhat red- 
dened objects (see the bottom right panel of Fig.|2j. For NGC 4594, 
the GCs are clearly discernable on the colour-magnitude diagram 
(CMD; see below), so we were able to set a more stringent colour 
cut: 0.5 < V—I < 1.5. We also applied a brightness cut, eliminat- 
ing GC candidates brighter than Gl, the most luminous M31 GC 
(Mv = -10.55, Richetal. 1996). 

To avoid complications related to excessive reddening, we ex- 
cluded candidates lying within the dust lanes of their parent galax- 
ies. Since this normally accounted for only a small fraction (~ 5%) 
of the total detected cluster population, this should have minimal 
effects on our global conclusions. However, we have attempted to 
correct this bias in the same manner as for our incomplete spatial 
coverage of the target galaxies (see below). 

V versus V — I CMDs for the candidate GCs in five of our 
program galaxies^ are shown in Fig.|2| together with that of the GC 
system of the Milky Way as a comparison. One rather unusual find- 
ing in this context is that the CMD of NGC 45 17 reveals a relatively 
large number of GC candidates with 0.3 < V — I < 0.6, which is 

2 For the CMDs and the photometry tables of the GCs in NGC 4565 and 
NGC 4594, we refer to Kissler-Patig et al. (1999) and Larsen, Forbes & 
Brodie (2001), respectively 
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Figure 2. V vs. V—I colour-magnitude diagrams of candidate GCs in the five sample galaxies for which GC photometry was not published before, compared 
with that of the GCs in the Milky Way as taken from the 2000 update of the compilation of Harris (1996; bottom right panel). All panels cover the same range 
in absolute V-band magnitudes under the assumption of the galaxy distances in Table fTI Filled symbols represent observed data, whereas open symbols in 
the bottom right panel represent dereddened data of the Milky Way GCs. The large reddening corrections for many Milky Way GCs prompted us to use fairly 
generous colour cuts in selecting GC candidates for spirals with HST pointings encompassing the dusty discs (see Section lZSl . 



bluer than Galactic GCs. However, the sizes of all these candidates 
(which are located on the WF chips of WFPC2) are fully consistent 
with GCs at the distance of NGC 45 17. Hence we choose to retain 
those sources as GC candidates. Spectroscopy will be needed to 
confirm the nature of these sources. 

The photometry and positions of the 50 brightest GC candi- 
dates in the five galaxies for which GC photometry was not pub- 
lished before are listed in Tables lBTllB5l (in appendixlHi. 



1981). The predicted star counts at the Galactic latitude and lon- 
gitude of each of the sample galaxies were subjected to the same 
colour and brightness cuts as the GC candidates and then subtracted 
from their respective bins. The total number of contaminating stars 
passing these criteria was small, generally ~ 3 — 5 stars per point- 
ing. 



3 RESULTS: PROPERTIES OF THE GLOBULAR 
CLUSTER SYSTEMS 

3.1 The total number of globular clusters 

The method we used to estimate the total number of GCs in the 
target galaxies was to make corrections for (i) contamination by 
Galactic foreground stars and (ii) spatial coverage. The latter was 
done by comparing the target GC systems to that of the Milky Way, 
as detailed below. 

3.1.1 Correction for Galactic foreground stars 

First, our estimates were corrected for foreground contamination, 
which was estimated from Galactic models (Bahcall & Soneira 



3.1.2 Correction for completeness 

Artificial star experiments were performed by using the DAOPHOT- 
II task ADDSTAR to estimate the completeness of the finding algo- 
rithm. These were carried out in the usual fashion, using artifical 
GCs generated from the PSFs constructed during the photometry 
measurements. These 'fake' GCs were added in groups of 100 in 
randomly placed positions on each chip. For the WF chips, the 50% 
completeness level was generally 1 .0 - 1 .5 mag beyond the turnover 
of the GC luminosity function (GCLF), ensuring accurate determi- 
nation of the turnover point and dispersion of the GCLF. The WF- 
chip completeness functions are shown in Fig.|3| For each magni- 
tude bin and each galaxy, a completeness fraction was calculated, 
and the number of objects in the bin was divided by this number to 
produce a completeness-corrected value. 
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Figure 3. V-band completeness functions for the WFPC2 photometry 
of GC candidates. From left to right (at the 50% completeness level), 
the curves represent the WF-chip completeness functions for NGC 7814, 
NGC 4594, NGC 4013, NGC 3628, NGC 4517. and IC 5176, respectively. 



from the McMaster catalog (Harris 1996), containing 141 GCs 
with appropriate {V — I) colour information. However, we adopted 
Nmw = 160 ± 20 (van den Bergh 1999) as the total number of 
GCs in the Milky Way. The undetected Milky Way GCs are as- 
sumed to be behind the Galactic bulge. Since the scale lengths and 
scale heights of the thin discs of the sample galaxies are very simi- 
lar to those of the Milky Way (e.g., Bothun, Harris & Hesser 1992; 
Morrison, Boroson & Harding 1994; Baggett, Baggett & Anderson 
1998; Gavazzi et al. 2000; Wu et al. 2002), the counterparts of such 
GCs would not be detected in our observations either. We there- 
fore assume that the percentage of obscured GCs is the same in all 
sample galaxies. 

We calculated a statistical uncertainty of the total number of 
GC candidates calculated this way by considering Poisson errors in 
the observed number of GCs, Poisson errors in the average number 
of Milky Way GCs in the mask (averaged over the four different 
orientations), and errors in the number of contaminating foreground 
objects. Our final numbers for the total population of GC candidates 
are listed in Table|3| 



3.1.3 Extrapolation over the luminosity function 



In preparation for the calculation of specific GC frequencies as de- 
fined by Harris &. van den Bergh (1981), the total number of GCs 
is defined as twice the number of GCs brighter than the turnover 
magnitude of the GCLF, where the GCLF is assumed to be a Gaus- 
sian (in magnitude units). We therefore fitted the completeness- 
corrected GCLFs by a Gaussian with a standard deviation (7 = 1.1 
mag (Harris 1996; Whitmore 1997) and derived the total number 
of clusters from the Gaussian fit. 



3.1.4 Correction for incomplete spatial coverage 

To estimate the total number of GCs in each galaxy, we need to 
correct for the incomplete spatial coverage of our data. Since each 
pointing only covered part of the galaxy, an extrapolation over the 
entire galaxy was needed to estimate the total GC population. In the 
case of elliptical and SO galaxies, this extrapolation is often done by 
evaluating the radial distribution of GCs within the observed radial 
extent of the galaxy. However, the low number statistics of GCs 
in most spiral galaxies does not allow this to be done accurately. 
Hence the correction for spatial coverage was done by means of a 
direct comparison to the positions of GCs in the Milky Way. Since 
this process has been described and illustrated in detail in Kissler- 
Patig et al. (1999), we will only summarize the process here. 

We created a mask defined by our spatial coverage of each 
galaxy, and applied this mask to the Milky Way. By calculating 
the number of objects we would detect in the Milky Way mask, if 
placed at the distance of the target galaxy, we calculated the total 
number of GCs Ncc as follows: 

M AT Nmw 
J-^GC ~ Nobs -r, 

l^MSK 

where Nmsk is the number of objects detected in the mask, and 
Nmw is the total number of GCs in the Milky Way. It should be 
noted that there are four different orientations of the mask that 
would preserve the position angle of the target galaxy (these can 
be seen as reflections of the mask across its horizontal or vertical 
axes). Our final value for Nmsk was the average of the values of 
all four possible orientations, and the standard deviation of the four 
values was taken into account in the uncertainty of Nmsk. 

For this study, data regarding Milky Way GCs were taken 



3.1.5 Comparison with previous studies 

A potential caveat of our method to correct for incomplete spatial 
coverage is the implicit assumption that the spatial distribution of 
the GCs in the sample galaxies is similar to that of the Milky Way 
globular cluster system. Hence a comparison with other methods 
is in order. Unfortunately, there have been only very few ground- 
based, wide-field studies of spiral galaxy GC systems, which is 
likely due to the complications discussed in the Introduction. Al- 
though this is slowly changing with the use of wide-field cameras 
on 4-m-class telescopes (Rhode & Zepf 2002), we can currently 
only make a useful comparison for one galaxy: NGC 7814, for 
which a previous (ground-based) study of the GC system was per- 
formed by Bothun et al. (1992). Bothun et al. found a much larger 
number of GC candidates than we did (Nqc = 498 ± 164 ver- 
sus our 106 ± 28). From an inspection of our WFPC2 frames, we 
suspect that the different result is due to their counting method, 
which was based on object overdensities (relative to a background 
field) rather than colours and sizes. We find a large number of small 
galaxies in our WFPC2 field, which Bothun et al. (1992) likely 
counted along with 'real' GCs, while they were excluded from our 
GC candidate list. This once again illustrates HSTs unique power 
in discriminating GCs from small, nucleated galaxies. 

A perhaps more useful comparison of our method with oth- 
ers can be done for the case of NGC 4594, for which Larsen et al. 
(2001) used the same HST data, but applied a different method: 
Larsen et al. constructed radial distribution functions of GCs to 
correct for incomplete spatial coverage. While that method might 
in principle be expected to produce more accurate results than the 
Milky Way mask method used here, it requires a large population 
of GCs to avoid complications related to small number statistics. 
(In our galaxy sample, this requirement is only met in the case 
of NGC 4594). In any case, a comparison of Larsen et al.'s re- 
sults with ours is very encouraging: They find a total population 
of iVoc = 1150 ±575, compared to iVcc = 1270 ±308 (statisti- 
cal error only) in this study. Especially heartening in this respect is 
that NGC 4594 (the Sombrero galaxy) is a typical Sa galaxy, which 
might perhaps not be expected to have a GC system particularly 
similar to that of the Milky Way, an Sbc galaxy. 
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3.2 Cluster specific frequencies 

Traditionally, the specific frequency of GCs in galaxies, Sn, is de- 
fined as Sjv = Ngc X ioO-4(«v'+i5) (Harris & van den Bergh 
1981), i.e., the total number of GCs per unit galaxy l/-band lumi- 
nosity (normalized to My = —15). Sn was introduced primar- 
ily for use in early-type (E and SO) galaxies, where there is little 
variation in stellar populations between individual galaxies. This 
is not necessarily the case for spirals, and indeed Sn is generally 
known to increase along the Hubble sequence (going from late- 
to early-type). Averaging specific frequencies per Hubble type, the 
recent compilation of Ashman & Zepf (1998) supplemented by the 
study of two spiral galaxies by Kissler-Patig et al. (1999) yields 
(Sn) = 0.4 ± 0.2 for Sc spirals to 1.9 ± 0.5 for elliptical galax- 
ies outside galaxy clusters (see Section|4}. Zepf & Ashman (1993) 
attempted to account for differences in stellar mass-to-light ratios 
among galaxy Hubble types in a statistical sense by introducing a 
parameter T (hereafter TzAga) to be the number of GCs per unit 
stellar mass (10^ M©) of a galaxy. Conversion from luminosity to 
mass was achieved my assuming a characteristic M / Lv value for 
each galaxy Hubble type. 

3.2.1 Total specific frequencies 

The ^jv and Tzaos values of the target galaxies are listed in Ta- 
ble|3| calculated using the absolute magnitudes given in TableQand 
the number of GCs derived in Section lTTI The quoted uncertainties 
were derived by taking into account the random errors discussed 
above as well as an uncertainty of 0.2 mag in the (total) absolute, 
dereddened V-band magnitudes of the galaxies. Tzaos values were 
calculated by converting the galaxy luminosities to masses follow- 
ing Zepf & Ashman (1993), using MjLy values of 5.4, 6.1, and 
5.0 for Sa, Sab-Sb, and Sbc-Sc galaxies, respectively (cf. Faber & 
Gallagher 1979). We note that the uncertainty in the distance of a 
galaxy has a negligible effect on the 5'jv and Tza93 values, since 
the change in the total number of GCs (due to a change in the physi- 
cal area covered by the WFPC2 frames) is compensated by a similar 
change of the luminosity of the galaxy. 

3.2.2 Bulge specific frequencies 

As already mentioned in the Introduction, a view is emerging that 
inner metal-rich GCs in spiral galaxies may be associated with 
their bulges rather than with their (thick) discs. Ample evidence 
for this is available for the Milky Way and M 3 1, where the metal- 
rich GC systems show kinematics, metallicities, and spatial dis- 
tributions matching those of the underlying bulge stars (Minniti 
1995; Cote 1999; Barmby et al. 2001; Perrett et al. 2002). Such 
a physical association is also likely present for elliptical galaxies, 
where the spatial distribution of metal-rich GCs typically follows 
that of the spheroidal galaxy light distribution very closely, whereas 
the metal-poor GC system typically has a more extended distribu- 
tion (e.g.. Ashman & Zepf 1998 and references therein). If indeed 
this association can be confirmed for a large sample of spirals as 
well, it would provide an important causal link between the for- 
mation of a spheroidal stellar system and that of metal-rich GCs. 
A suggestion that this might be the case was recently provided by 
Forbes et al. (2001) who compared the GC systems of the Milky 
Way (an Sbc galaxy), M 3 1 (Sb), and M 104 (Sa). They argued that 
the "bulge specific frequency (bulge Sn)" - which they defined as 
the number of metal-rich GCs within 2 bulge effective radii from 
the galaxy centres divided by the bulge luminosity (normalized to 



Mv = —15) - was consistent among those three galaxies (bulge 
5'jv ~ 0.5), and similar to typical values for field elliptical galaxies 
as well. 

Structural bulge properties and bulge-to-total luminosity ra- 
tios (hereafter B /T ratios) for the sample galaxies are listed in 
TableQ For NGC 4517, NGC 4565, NGC 4594 and NGC 7814, 
these numbers were taken from the literature (Gavazzi et al. 2000; 
Wu et al. 2002; Baggett, Baggett & Anderson 1998; Bothun, Har- 
ris & Hesser 1992, respectively). The bulge/disc decompositions 
for the other sample galaxies are described in Appendix |^ For 
NGC 4013 and IC 5176, we used our own WFPC2 images (in the 
F555W band). For NGC 3628, we used the WIYN V-band im- 
age published by Howk & Savage (1999) which was graciously 
made available to us. Among the sample galaxies, only the bulges 
of NGC 4013, NGC 4517, NGC 4594 and NGC 7814 were well fit 
by a de Vaucouleurs' profile. The bulges of NGC 3628, NGC 4565 
and IC 5176 were much better fit by an exponential profile. Half- 
light radii ri/2 for the bulges of the latter galaxies were calculated 
as 1.678 ro, where ro is the exponential scale length of the profile. 
The final half-light bulge radii ri/2 of all sample galaxies are listed 
in the last row of TableQ 

To identify 'metal-rich' GCs in our sample galaxies, we con- 
sidered a lower limit of [Fe/H] = —1, which is the location of 
the 'dip' in the metallicity distribution of the Milky Way GC sys- 
tem (Harris 1996). This is shown in Fig.|4] superposed onto the 
GC metallicity histograms of the sample galaxies. The metallic- 
ities were derived by transforming the V — I colours (corrected 
for Galactic extinction) into [Fe/H] values using the recent cali- 
bration of Kissler-Patig et al. (1998). The B- 7 colours of GCs in 
NGC 4565 were converted into [Fe/H] values using the relation for 
Milky Way GCs, given in Couture, Harris & Allwright (1990). Note 
that the [Fe/H] values for GC candidates in external spiral galaxies 
as derived from colours are likely overestimates given the possibil- 
ity of reddening by dust, especially in the inner regions. Thus, the 
number of "metal-rich" clusters resulting from this exercise has to 
be regarded as an upper limit. 

To allow a direct comparison with the 'bulge 5'jv' values de- 
fined in Forbes et al. (2001), we counted the total number of 'metal- 
rich' GC candidates as well as those with a projected galactocentric 
radius < 2 ri/2 identified in our images. Correction for incomplete 
spatial coverage was done by dividing the observed number of GCs 
by the fraction of the total area within 2ri/2 sampled by our im- 
ages. Since the surface density of GCs typically falls off with in- 
creasing galactocentric radius, this correction procedure was done 
in elliptical annuli, i.e., multiplying the number of GCs in each an- 
nulus by the ratio of total annulus area to that which was actu- 
ally observed. We chose to skip this measurement for the case of 
NGC 3628 since the HST observations of NGC 3628 only covered 
a very small portion of the inner bulge, rendering the extrapolated 
result very uncertain. In the cases of NGC 4565 and IC 5176, we 
only counted GCs on one side of the disc for this purpose, namely 
the side that was visually not significantly impacted by dust ex- 
tinction. The bulge 5jv values for the sample galaxies were then 
calculated using the B /T ratios listed in Table Q The results are 
listed in Table|3] 

The bulge 5jv values within 2ri/2 for most of our galaxies 
are roughly consistent with a value of ~0. 1, with the exception 
of NGC 4594, the Sombrero galaxy, for which we find a value 
of ~0.5. At face value, the low bulge Sjv values for most of our 
galaxies seem to be inconsistent with the prediction of Forbes et 
al. (2001) who suggested that bulge 5jv values within 2r-i/2 are 
constant among spiral galaxies with a value around 0.5. However, 
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Table 3. Number of GCs and specific frequencies for the sample galaxies. 



Galaxy 


GCs 


GCs 


Total 5jv 


TzAOS 


Metal-rich GCs 


Bulge Sjv 


Metal-rich GCs 


Bulge Sm 




(detected) 


(total) 






(all radii) 


(all radii) 


(r < 2ri/2) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



NGC 3628 


92 ± 10 


497 ± 110 


1.9 ±0.2 


3.7 ± 1.0 


213 ±48 


0.83 ±0.11 


a 




NGC 4013 


69 ±8 


243 ± 51 


1.1 ±0.3 


2.2 ± 0.7 


95 ±20 


0.44 ±0.18 


3 ± 3 


0.06 ± 0.06 


NGC 4517 


62 ±8 


270 ± 60 


0.6 ± 0.2 


1.4 ±0.5 


81 ± 18 


0.18 ±0.16 


b 




NGC 4565 


40 ±6 


204 ± 38= 


0.6 ± 0.2 


1.0 ±0.3 


122 ± 23 


0.31 ±0.15 


20 ±5'' 


0.17 ±0.15 


NGC 4594 


159 ± 13 


1270 ± 308 


1.7 ±0.6 


3.6 ± 1.1 


691 ± 167 


0.91 ±0.27 


429 ± 214 


0.77 ± 0.44 


IC 5176 


57 ± 8 


132 ± 25 


0.5 ±0.1 


1.1 ±0.3 


67 ± 12 


0.25 ± 0.09 


7 _l_ 2d 


0.07 ± 0.07 


NGC 7814 


17 ±4 


106 ± 28 


0.7 ± 0.2 


1.5 ±0.5 


31 ±9 


0.20 ± 0.09 


8 ± 5 


0.06 ± 0.06 



Notes. Column (1): Galaxy name. Column (2): Detected number of GC candidates in WFPC2 images of galaxy. Column (3): Total number of GC 
candidates around galaxy. Column (4): Total specific frequency of GC candidates. Column (5): T-parameter as defined by Zepf & Ashman (1993) as the 
number of GCs per unit galaxy mass (in terms of 10® Mq). Column (6): Total number of metal-rich GCs ([Fe/H] > —1) around galaxy. Column (7): 
Total bulge specific frequency of metal-rich GC candidates. Column (8): Number of metal-rich GC candidates within 2 bulge half-light radii. Column 
(9): Bulge specific frequency of metal-rich GC candidates within 2 bulge half-Hght radii. Note a: Not measured due to too small spatial coverage of inner 
bulge. Note b: Bulge radius too small to find bulge GCs. Note c: Taken from Kissler-Patig et al. (1999). Note d: Number is twice the GC candidates found 
on the non-dusty side of the nucleus. 
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Figure 4. Metallicity distributions of GC candidates in the spiral galaxies in our sample, and of GCs in the Milky Way for comparison. The y-axis scaling is 
different for each galaxy, but the small tickmarks always represent 1 GC. Open histograms represent the total observed GC system, and hatched histograms 
represent observed GCs within 2 bulge half-light radii from the galactic center. The vertical dashed line at [Fe/H] = — 1 separates 'metal-poor' from 'metal-rich' 
GC candidates. 



systematic uncertainties of the inner bulge Sn values as derived 
from our data can be substantial. E.g., we generally optimized our 
observing strategy to detect GCs in the halos of the target galaxies 
(cf. Fig.0, and dust extinction effects are strong in the innermost 
bulge regions. Hence we defer further discussion on the universal- 
ity of the bulge Sjv until high-resolution near-IR data are available. 
Our following discussion on trends of Sjv values along the Hubble 
sequence will therefore focus on the 'total' Sjv values. 



4 DISCUSSION 

In the context of the 'major merger' (Ashman & Zepf 1992) or 
'multi-phase collapse' (Forbes, Brodie & Grillmair 1997) scenar- 
ios for the formation of early-type galaxies, the higher Sn values 
of ellipticals (with respect to spirals) are due to 'extra' GCs that 
were formed with high efficiency during the event that also formed 
the spherical stellar component (bulge). If this is indeed how bulges 
of spirals form in general, one would expect Sn to vary systemat- 
ically along the Hubble sequence of spirals (increasing from late- 
type to early-type spirals). If instead bulges were formed from disc 
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stars by means of redistribution of angular momentum through mi- 
nor perturbations and/or bar-like instabilities as in the scenario of 
Pfenniger & Norman (1990), one would not expect to see any sig- 
nificant variation along the Hubble sequence, since no GCs would 
be formed in this secular building of bulges. 

Since our galaxy sample spans a large range of Hubble types 
(from Sa to Sc), we are in a good position to test this. We plot 
the Sn and Tzaos values versus Hubble type and B /T ratio in 
Figure |5| Values for the Milky Way and M31 (taken from Forbes 
et al. 2001 and references therein) are also included in that Figure. 
There are a number of remarkable aspects of these results, which 
we discuss in order below. 

First of all, there does not seem to be any significant difference 
between trends derived from Tza9s and those from Sn- Hence, we 
will only discuss trends involving Sn, with the understanding that 
the trends involving Tza9s are consistent with them. 

A general trend in Figure |5| is that Sn stays basically con- 
stant at a value of ~0.55 for spirals with B/T ^ 0.3 (roughly 
corresponding to Hubble types Sb and later). As these galaxies are 
clearly dominated in mass by their halos, this constant Sn (and 
TzAgs) suggests that the formation process of massive (late-type) 
spiral galaxies involves the creation of a quite constant number of 
GCs per unit galaxy luminosity (and mass). We suggest that this 
population of GCs represents a 'universal', old halo population that 
is present around each galaxy. 

It is now well established that elliptical galaxies in environ- 
ments similar to that of the spirals in our sample (i.e., outside 
galaxy clusters) have higher specific frequencies than these late- 
type spirals. A x^-weighted average Sn value for ellipticals out- 
side galaxy clusters published in the literature is Sjv = 1.9 ± 0.6 
(using Sjv data from Ashman & Zepf 1998; Brown et al. 2000; 
Georgakakis, Forbes & Brodie 2001; Goudfrooij et al. 2001). If in- 
deed the higher values of Sn for elliptical galaxies are due to an 
'extra' population of GCs which was formed during the formation 
process of spheroidal stellar systems (such as bulges) involving a 
rapid dissipative collapse in which the physical conditions are such 
that giant star clusters are formed very efficiently (i.e., as in the 
'merger' and 'multi-phase collapse' scenarios mentioned above), 
then one would naively expect the Sn values to scale linearly with 
B/r ratio from 5]v = 0.55±0.25 at B/T = OtoS'jv = 1.9±0.6 
at B /T — 1. This area is hatched in the middle right panel of Fig.|5| 

Interestingly, the Sn values of most sample galaxies are in- 
deed consistent with this simple prediction. There are however two 
galaxies for which this is clearly not the case: NGC 3628 (Sb) and 
NGC 78 14 (Sab). 

The result for NGC 3628 is quite remarkable. This Sb galaxy 
turns out to have a Sn value that is significantly higher than all 
other Sb or She galaxies with Sn measurements to date. How- 
ever, NGC 3628 is a peculiar galaxy. It is part of the Leo Triplet 
(Arp 317; Arp 1966). Large-area optical observations reveal a tidal 
plume extending ~ 100 kpc towards the east, as well as a bridge 
between NGC 3628 and NGC 3627 (Burkhead & Hutter 1981; 
Chromey et al. 1998), both of which contain large amounts of Hi 
gas (5.4x10* Mq, which is 15% of the total H I mass observed 
in the main body of NGC 3628; Haynes, Giovanelli & Roberts 
1979). Rots (1978) constructed three-body orbital models for the 
tidal interactions between NGC 3627 and NGC 3628 which repro- 
duce the formation of the plume and bridge, yielding an age for 
the plume of 8 x 10* yr since perigalacticon. This age was con- 
firmed by photometric measurements of bright clumps in the plume 
by Chromey et al. (1998). It is therefore quite conceivable that the 
high Sn value of NGC 3628 may be due to an 'extra' supply of 



star clusters which were formed during the gaseous interaction that 
also caused the tidal plume and bridge. This situation seems to have 
occurred in M82, another nearby starburst galaxy which has under- 
gone a past starburst triggered by tidal interactions with a gas-rich 
neighbour, during which a population of GCs did indeed form (see, 
e.g., de Grijs, O'Connell & Gallagher 2001). As the gas from which 
such second-generation star clusters are formed is expected to be of 
relatively high metallicity (relative to halo GCs), this picture is con- 
sistent with the observation that the (total) bulge Sn of NGC 3628 
(i.e., the number of metal-rich GC candidates per unit bulge lu- 
minosity) is higher than those of the other Sb-Sc galaxies in our 
sample by a factor similar to the total Sn- This idea should be ver- 
ified by performing wide-field imaging as well as spectroscopy of 
a significant number of GC candidates in NGC 3628. 

As to the deviant result for NGC 7814, it is perhaps somewhat 
surprising to find that the two earliest-type spirals in this sample 
(NGC 4594 and NGC 7814) have very different specific frequen- 
cies. NGC 4594 (Sa) has Sn values that are similar to those of field 
ellipticals (e.g.. Ashman & Zepf 1998) and SO galaxies (Kundu & 
Whitmore 2001), and thus seems to have formed GCs at a similarly 
high efficiency relative to its field stars. It is therefore highly un- 
likely that bulges like the one of NGC 4594 can be formed through 
the 'secular evolution' scenario of Pfenniger & Norman (1990). On 
the other hand, NGC 7814 has a much lower Sn, which is consis- 
tent with the values found for the bulk of the latest-type spirals. 
This difference in 5'jv persists when considering the 'bulge Sn' 
values. At face value, this result suggests that the (star-) forma- 
tion histories of galaxy bulges of early-type spirals can be signifi- 
cantly different from one galaxy to another. The formation process 
of some bulges (such as that of NGC 4594) may have involved 
dissipative collapse, being accompanied by the creation of a sig- 
nificant number of (metal-rich) GCs, while bulges in other spiral 
galaxies (like NGC 7814) may have been built from the redistribu- 
tion of angular momentum through bar-like instabilities as in the 
'secular evolution' scenario. In this context, it may well be relevant 
that NGC 7814 is the least luminous (and thus presumably least 
massive) spiral galaxy in our sample, being ~ 5 times less lumi- 
nous than NGC 4594. As Pfenniger (1991) showed, the accretion 
of small satellite galaxies (with a mass of ~ 10% of that of the mas- 
sive galaxy) can trigger the formation of a bar and produce a bulge 
through secular evolution. It is tempting to suggest that the differ- 
ence in Sn value between NGC 4594 and NGC 7814 is simply 
due to the fact that satellite galaxies at a satellite-to-giant galaxy 
mass ratio of ~10% are more common at low (giant) galaxy lu- 
minosities than at high luminosities. Thus suggestion is certainly 
consistent with observed luminosity functions of galaxies in poor 
groups (Jerjen, Binggeli & Freeman 2000; Flint, Bolte & Mendes 
deOliveira 2001). 

Given the importance of this issue in the context of our un- 
derstanding of the formation of the Hubble sequence of galaxies, 
we argue that this suggestion be tested with a statistically signifi- 
cant sample of (edge-on) Sa-Sab spirals which covers an appro- 
priate range of luminosities. This can now be realized with only 
very modest amounts of observing time using the recent ACS cam- 
era aboard HST, which reaches some 2 magnitudes deeper than 
WFPC2 at a given exposure time. The current generation of large- 
area CCD mosaics on ground-based 4-m-class telescopes will be 
very valuable as well in terms of covering the whole GC popula- 
tions around such galaxies in one shot. 
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Figure 5. Relationships of the specific frequency parameters (^jv and bulge ^jv) and the T parameter (introduced by Zepf & Ashman 1993; hence we identify 
it as TzA9s) with galaxy Hubble type and bulge-to-total luminosity ratio for the sample galaxies as well as the Milky Way and M31 (data for the latter two 
galaxies were taken from Forbes et al. 2001). The identifications of the galaxies are plotted next to the accompanying data points (at a location least impacting 
the other symbols). The hatched area in panel (d) represents the region of expected 5]v values for a simple model described in Section|3 In panels (e) and (f) 
(i.e., the bulge Sn plots), the filled circles represent the 'total' bulge ^jv values, while the open circles represent the bulge ^jv values within 2 half-light radii 
of the bulges (see Section lX2l 2). To avoid overcrowding, the galaxy identifications are only plotted along with the filled circles. 



5 SUMMARY 

We have studied the GC systems of 7 giant, edge-on spiral galaxies 
covering the Hubble types Sa to Sc, allowing us to study the vari- 
ation of the properties of GC systems along the Hubble sequence. 
Our high-resolution HST/WFPC2 data (supplemented by archival 
WFPC2 data for two galaxies) reached ~ 1.5 mag beyond the turn- 
over magnitude of the GC luminosity function for each galaxy. 

Specific frequencies of GCs {Sn values) were evaluated by 
comparing the numbers of GCs found in our WFPC2 pointings with 



the number of GCs in our Milky Way which would be detected 
in the same spatial region if placed at the distance of the target 
galaxies. Results from this method were compared with the more 
commonly used method of constructing radial distribution func- 
tions of GCs to correct for incomplete spatial coverage in the case 
of NGC 4594 (where both methods were possible), and found to be 
consistent with one another to within la. 

The Sn values of spirals with B/T ^ 0.3 (i.e., spirals with 
a Hubble type later than about Sb) are consistent with a value of 
Sn ~ 0.55 ± 0.25. We suggest that this population of GCs repre- 
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sents a 'universal', old halo population that is present around each 
galaxy. 

Most galaxies in our sample have Sn values that are consis- 
tent with a scenario in which GC systems are made up of ( i) the 
aforementioned halo population plus (ii) a population that is associ- 
ated with bulges of spirals, which grows linearly with the luminos- 
ity (and mass) of the bulge. Such scenarios include the 'merger sce- 
nario' for the formation of elliptical galaxies as well as the 'multi- 
phase collapse' scenario, but it seems inconsistent with the 'secular 
evolution' scenario of Pfenniger & Norman (1990), in which bulges 
are formed from disc stars by means of the redistribution of angular 
momentum through minor perturbations and/or bar instabilities. 

On the other hand, the bulge-dominated spiral NGC 7814 
shows a low Sn value, consistent with those of the latest-type 
spirals. NGC 7814 is the least luminous galaxy in our sample. 
Based on observed luminosity functions of galaxies in poor groups, 
we suggest that the 'secular evolution' scenario to build bulges in 
early-type spirals is most viable for low-luminosity spirals. 

Thus, our results suggest that the formation histories of galaxy 
bulges of early-type spirals can be significantly different from one 
galaxy to another. Given the importance of our understanding of 
the formation of the Hubble sequence of galaxies in the context of 
galaxy evolution, we argue that the GC systems of a statistically 
significant sample of luminous, edge-on Sa- Sab spirals be studied 
in the near future. 
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APPENDIX A: BULGE/DISC DECOMPOSITIONS 

Since bulge/disc decompositions are not the main topic of this pa- 
per, we will only provide a brief description of our method here. 
A more detailed description will be provided elsewhere. We first 
determined the scale heights of the discs at 4 independent positions 
along their major axes, outside the region where the bulge contribu- 
tion dominates. In general, the positions of the galaxy planes were 
determined by assuming symmetrical light distributions with re- 
spect to the disc plane, and the vertical profiles on either side of the 
plane were averaged together. However, the latter was not done for 
the bulge fits on NGC 4013 and IC 5176, since the WFPC2 frames 
only covered a large enough spatial extent on one side of the galaxy 
centre. 

The vertical brightness distributions of the discs were well fit- 
ted by an exponential profile 

Ed(2:) = ED,oexp(-z/2o), (Al) 

where zo is the disc scale height. We found that as long as one 
stays outside the central areas where there is an obvious contribu- 
tion from bulge light, the disc scale heights generally do not exhibit 
any significant variation with distance along the major axis. This is 
consistent with the findings of de Grijs, Peletier & van der Kruit 
(1997). The average disc scale heights are listed in TableQ 

To obtain structural parameters for the bulges of these galax- 
ies, we used a least-square program that fits a two-component 
(disc-l-bulge) model to the minor-axis brightness profiles. The disc 
component was fit by an exponential profile and the bulge com- 
ponent was fitted by two different functions, namely ( i) a de Vau- 
couleurs' (1953) profile 

Eb(z) = Eb,c acp{-imS{{z/z^f^* - 1)), (A2) 

where Zc is the effective radius along the minor axis and Eb.c is 
the intensity &t z = Zc, and (ii) a (second) exponential profile. For 
each galaxy, the disc scale height zq as found above was fixed in 
the fit. Areas where dust extinction is prominent (such as the cen- 
tral dust lanes) were flagged and ignored during the fitting process. 
The choice between the de Vaucouleurs' profile and the (second) 
exponential profile was dictated by the value of the fit. 

Among the sample galaxies, only the bulges of NGC 4013, 
NGC 4517, NGC 4594 and NGC 7814 were well fit by a de Vau- 
couleurs' profile. The bulges of the other galaxies were better fit 
by an exponential profile. The effective radii (or exponential scale 
radii rcxp) of the bulges are listed in TableQ after converting them 
to equivalent radii r = yfab — 6/\/l — where a and b are the 
semimajor and semiminor axes of the ellipse, and e its ellipticity. 
For NGC 3628, NGC 4013, NGC 4517 and IC 5176, the bulge el- 
lipticities were determined by using the ellipse fitting program in 
the ISOPHOTE package within IRAF. The dusty areas of the galax- 
ies were flagged and ignored in the ellipse fitting process. For the 
galaxies with exponential bulge profiles, we also list the bulge half- 
light radii (= 1.678 rexp) for comparison with the effective radii of 
the bulges with de Vaucouleurs' profiles. 



APPENDIX B: TABLES OF GLOBULAR CLUSTER 
CANDIDATES 

The photometry and astrometry of the 50 brightest GC candidates 
on the WFPC2 images of the galaxies for which GC photometry 
was not published before are given in Tables lBn - IB5l For the pho- 
tometry of GC candidates in NGC 4565 and NGC 4594 we refer to 
Kissler-Patig et al. (1999) and Larsen et al. (2001), respectively. 



Table Bl. Photometry and astrometry of the 50 brightest globular cluster 
candidates on the WFPC2 frames of NGC 3628. The object list is sorted on 
V magnitude (brightest first). 



ID 


RA 
(J2000) 


DEC 
(J2000) 




V 
mag 






V-I 
mag 




3628-1 


11 


20:14.64 


+13 


32:27.6 


20 


65 ±0 


01 





41 ±0 


02 


3628-2 


11 


20:21.12 


+13 


35:38.4 


20 


97 ±0 


02 





92 ±0 


03 


3628-3 


1 1 


20: 19.92 


+13 


35:45.6 


21 


83 ±0 


02 


1 


01 ±0 


03 


3628^ 


1 1 


20-19 44 


+13 


35-50 4 


21 


95 ±0 


01 





93 ±0 


03 


3628-5 


1 1 


20' 14 40 


+13 


33-57 6 


21 


96 ±0 


03 





87 ±0 


04 


3628-6 


11 


20:17.76 


+13 


35:45.6 


22 


07 ±0 


03 


1 


21 ±0 


04 


3628-7 


11 


20:17.76 


+13 


37:22.8 


22 


10 ±0 


02 





81 ±0 


03 


3628-8 


1 1 


20:12.72 


+13 


34:04.8 


22 


12 ±0 


24 





44 ±0 


24 


3628-9 


I \ 


90-1 S (lO 


+13 


33-32 4 


22 


15 ±0 


02 


1 


01 ±0 


04 


3628-10 


1 1 




+13 


33-36 


22 


23 ±0 


02 





87 ±0 


06 


3628-1 1 


11 


20:12.72 


+13 


33:57.6 


22 


29 ±0 


08 





90 ±0 


11 


3628-12 


11 


20:19.92 


+13 


35:31.2 


22 


37 ±0 


05 





76 ±0 


08 


3628-13 


1 1 


nrt. 1 /I OA 


+13 


35:34.8 


22 


37 ±0 


05 


1 


49 ±0 


05 


3628-14 




20' 16 32 


+13 


35'38 4 


22 


44 ±0 


03 


1 


27 ±0 


03 


3628-15 




90- 1 R no 


+13 


35-31 2 


22 


50 ±0 


07 





87 ±0 


08 


3628-16 


11 


20:16.80 


+13 


35:52.8 


22 


53 ±0 


04 


1 


36 ±0 


04 


3628-17 


11 


20:17.04 


+13 


35:31.2 


22 


71 ±0 


06 


1 


66 ±0 


06 


3628-18 


1 1 


20: 15.84 


+13 


35:38.4 


22 


71 ±0 


04 


1 


29 ±0 


05 


3628-19 




20' 15 84 


+13 


35-34 8 


22 


72 ±0 


04 





97 ±0 


07 


3628-20 


I \ 


20-18 48 


+13 


36-07 2 


22 


73 ±0 


02 


1 


16 ±0 


03 


3628-21 


11 


20:17.76 


+13 


35:45.6 


22 


80 ±0 


04 





94 ±0 


06 


3628-22 


11 


20:16.08 


+13 


35:34.8 


22 


86 ±0 


10 





69 ±0 


13 


3628-23 


1 1 


20. 1 /.2o 


+13 


35:42.0 


22 


90 ±0 


04 


1 


08 ±0 


05 


3628-24 


I \ 


20-13 44 


+13 


34-04 8 


22 


96 ±0 


17 





66 ±0 


18 


3628-25 


I \ 


90-1 ^^ OS 

Z.\J. i u.uo 


+13 


35-38 4 


22 


99 ±0 


05 





73 ±0 


07 


3628-26 


11 


20:20.40 


+13 


34:40.8 


23 


04 ±0 


05 





70 ±0 


09 


3628-27 


11 


20:17.28 


+13 


34:44.4 


23 


06 ±0 


11 





94 ±0 


13 


3628-28 


1 1 


20:1 0.24 


+13 


35:27.6 


23 


08 ±0 


06 


1 


57 ±0 


06 


3628-29 


I \ 


90-1 (^ so 


+13 


36-36 


23 


11±0 


02 


1 


46 ±0 


04 


3628-30 


[ I 


20' 16 32 


+13 


35-31 2 


23 


11±0 


06 





68 ±0 


10 


3628-31 


11 


20:17.28 


+13 


35:38.4 


23 


12 ±0 


07 





66 ±0 


10 


3628-32 


11 


20:19.20 


+13 


35:34.8 


23 


15 ±0 


09 





81 ±0 


11 


3628-33 


1 1 


20:19.92 


+13 


35:52.8 


23 


18 ±0 


14 





70 ±0 


15 


3628-34 


[ I 


90- 1 R 00 


+13 


35-27 6 


23 


26 ±0 


10 


1 


78 ±0 


11 


3628-35 


1 1 


90-90 4.0 


+13 


35-34 8 


23 


26 ±0 


08 





46 ±0 


11 


3628-36 


11 


20:16.56 


+13 


35:42.0 


23 


30 ±0 


06 


1 


28 ±0 


07 


3628-37 


11 


20:17.04 


+13 


35:27.6 


23 


31 ±0 


08 


1 


37 ±0 


11 


3628-38 


11 


20:16.56 


+13 


35:34.8 


23 


37 ±0 


08 


1 


25 ±0 


10 


3628-39 


11 


20:20.40 


+13 


37:37.2 


23 


37 ±0 


03 


1 


86 ±0 


04 


3628^0 


11 


20:16.32 


+13 


35:31.2 


23 


39 ±0 


06 


1 


08 ±0 


09 


3628^1 


11 


20:18.96 


+13 


35:27.6 


23 


40 ±0 


17 





97 ±0 


19 


3628^2 


11 


20:18.72 


+13 


34:33.6 


23 


44 ±0 


03 


1 


63 ±0 


04 


3628^3 


11 


20:18.48 


+13 


33:43.2 


23 


45 ±0 


03 


1 


00 ±0 


04 


3628-44 


11 


20:17.04 


+13 


35:27.6 


23 


59 ±0 


10 


1 


93 ±0 


11 


3628-45 


11 


20:18.48 


+13 


36:25.2 


23 


61 ±0 


15 





90 ±0 


18 


3628-46 


11 


20:17.04 


+13 


35:31.2 


23 


62 ±0 


16 


1 


75 ±0 


17 


3628-47 


11 


20:17.52 


+13 


35:34.8 


23 


68 ±0 


08 





84 ±0 


11 


3628^8 


11 


20:14.88 


+13 


36:32.4 


23 


69 ±0 


06 





50 ±0 


12 


3628-49 


11 


20:15.84 


+13 


37:01.2 


23 


72 ±0 


17 





34 ±0 


21 


3628-50 


11 


20:16.80 


+13 


35:38.4 


23 


80 ±0 


10 


1 


44 ±0 


12 
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Table B2. Photometry and astrometry of the 50 brightest globular cluster 
candidates on the WFPC2 frames of NGC 4013. The object list is sorted on 
V magnitude (brightest first). 



ID 


RA 
(J2000) 


DEC 
(J2000) 




V 
mag 






V-I 
mag 




4013-1 


11 




443:56:24.0 


21 


30 ±0 


01 


1 


05 ±0 


04 


4013-2 


1 1 


CO O/l f\/' 

58:30.96 


+43:56:34.8 


21 


78 ±0 


02 


1 


41 ± 


03 


4013-3 


I \ 


58:29.52 


+43:56:02.4 


21 


88 ±0 


02 





91 ± 


03 


4013-4 


11 


58:33.36 


+43:57:18.0 


22 


12 ±0 


02 





80 ±0 


03 


4013-5 


11 


58:35.28 


+43:57:46.8 


22 


20 ±0 


01 





67 ±0 


02 


4013-6 


11 


CO 

58:29.76 


+43:56:52.8 


22 


25 ±0 


02 


1 


41 ±0 


03 


4013-7 


1 1 


58:32.64 


+43:56:16.8 


22 


36 ±0 


04 


1 


08 ± 


06 


4013-8 


1 1 


58-30 24 


+43:57:14.4 


22 


44 ±0 


02 


1 


00 ± 


03 


4013-9 


11 


58:31.44 


+43:56:13.2 


22 


45 ±0 


02 


1 


12 ± 


05 


4013-10 


11 


58:33.84 


+43:57:25.2 


22 


51 ±0 


03 





78 ±0 


05 


4013-11 


11 


58:28.56 


+43:56:52.8 


22 


55 ±0 


02 


1 


00 ±0 


04 


4013-12 


11 


58:31.68 


+43:57:00.0 


22 


69 ±0 


05 


1 


22 ±0 


07 


4013-13 


11 


5S-33 36 


+43:56:45.6 


22 


74 ±0 


02 


1 


10 ±0 


03 


4013-14 


11 


58:34.08 


+43:56:24.0 


22 


76 ±0 


10 





91 ± 


16 


4013-15 


11 


58:25.44 


+43:57:43.2 


22 


80 ±0 


10 





56 ±0 


11 


4013-16 


11 


58:31.44 


+43:58:08.4 


22 


83 ±0 


02 





59 ±0 


03 


4013-17 


11 


58:30.00 


+43:57:21.6 


22 


83 ±0 


01 





83 ±0 


03 


4013-18 


11 


58:29.52 


+43:56:56.4 


22 


88 ±0 


04 





96 ±0 


05 


4013-19 


11 


58:32.40 


+43:56:34.8 


22 


89 ±0 


04 


1 


08 ± 


05 


4013-20 


11 


58:27.36 


+43:56:09.6 


22 


98 ±0 


03 





94 ±0 


08 


4013-21 


11 


CO O o o 

58:32.88 


+43:56:27.6 


23 


02 ±0 


04 


1 


18 ±0 


06 


4013-22 


11 


58:30.96 


+43:57:00.0 


23 


06 ±0 


06 


1 


23 ±0 


07 


4013-23 


11 


58'33 36 


+43:55:51.6 


23 


18 ±0 


10 





95 ±0 


13 


4013-24 


11 


58:31.68 


+43:57:10.8 


23 


18 ±0 


04 





71 ±0 


06 


4013-25 


11 


58:29.76 


+43:57:10.8 


23 


29 ±0 


08 





85 ±0 


09 


4013-26 


11 


58:31.20 


+43:57:57.6 


23 


32 ±0 


03 





77 ±0 


05 


4013-27 


11 


CO O o o 

58:32.88 


+43:56:34.8 


23 


33 ±0 


05 


1 


09 ±0 


08 


4013-28 


11 


58'98 «n 

JO.^O.OU 


+43:56:13.2 


23 


34 ±0 


05 


1 


07 ±0 


07 


4013-29 


11 


58:32.40 


+43:57:57.6 


23 


44 ±0 


04 


1 


07 ±0 


04 


4013-30 


11 


58:31.20 


+43:57:03.6 


23 


46 ±0 


06 


1 


01 ±0 


08 


4013-31 


1 1 


58:32.64 


+43:56:42.0 


23 


54 ±0 


06 


1 


35 ± 


08 


4013-32 


11 


CO o 1 r*^ 

58:31.92 


+43:56:42.0 


23 


62 ±0 


06 


1 


36 ±0 


12 


4013-33 


11 


58:29.52 


+43:56:49.2 


23 


66 ±0 


06 


1 


05 ±0 


10 


4013-34 


11 


58:24.72 


+43:57:25.2 


23 


71 ±0 


03 





64 ±0 


04 


4013-35 


11 


58:30.72 


+43:57:39.6 


23 


77 ±0 


02 


1 


47 ±0 


06 


4013-36 


1 1 


CO O A 

58:33.84 


+43:57:32.4 


23 


77 ±0 


04 





91 ± 


05 


4013-37 


11 


CO -1 ^ f\ 

58:31.20 


+43:56:16.8 


23 


78 ±0 


04 





95 ±0 


07 


4013-38 


11 


58:30.72 


+43:56:34.8 


23 


81 ±0 


07 


1 


08 ±0 


12 


4013-39 


11 


58:33.84 


+43:56:20.4 


23 


96 ±0 


05 


1 


13 ±0 


09 


4013^0 


11 


58:28.56 


+43:57:03.6 


23 


97 ±0 


08 





87 ±0 


11 


4013^1 


11 


58:30.00 


+43:57:07.2 


23 


99 ±0 


05 





82 ±0 


09 


4013^2 


11 


58:32.40 


+43:57:14.4 


24 


04it0 


06 





76 ±0 


09 


4013^3 


11 


58:30.24 


+43:56:56.4 


24 


08 ±0 


09 


1 


08 ±0 


12 


4013^ 


11 


58:33.12 


+43:56:38.4 


24 


10 ±0 


08 


1 


10 ±0 


10 


4013^5 


11 


58:31.68 


+43:58:15.6 


24 


22 ±0 


10 


1 


04 ±0 


11 


4013-46 


11 


58:28.80 


+43:56:52.8 


24 


25 ±0 


08 


1 


02 ± 


12 


4013-47 


11 


58:34.08 


+43:56:42.0 


24 


28 ±0 


09 





88 ± 


12 


4013-48 


11 


58:33.84 


+43:56:45.6 


24 


33 ±0 


07 





98 ± 


13 


4013^9 


11 


58:31.92 


+43:56:34.8 


24 


36 ±0 


06 





99 ±0 


11 


4013-50 


11 


58:33.84 


+43:57:28.8 


24 


43 ±0 


05 





72 ±0 


10 
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Table B3. Photometry and astrometry of the 50 brightest globular cluster 
candidates on the WFPC2 frames of NGC 4517. The object list is sorted on 
V magnitude (brightest first). 



ID RA DEC V V-I 





(J2000) 


(J2000) 




mag 






mag 




4517-1 


12:32:43.68 


+00:06:07.2 


21 


70 ± 


02 





86 it 


02 


4517-2 


12:32:47.52 


+00:06:00.0 


21 


91 ± 


02 


1 


24 ± 


02 


4517-3 


12:32:43.68 


+00:06:32.4 


22 


05 ± 


03 


1 


20 ± 


03 


4517^ 


12:32:42.48 


+00:05:16.8 


22 


07 ±0 


05 





88 ± 


06 


4517-5 


12:32:48.72 


+00:05:31.2 


22 


10 ±0 


02 


1 


87it0 


03 


4517-6 


12-32-4fi 32 




22 


13 ± 


02 


1 


90 it 


11 


4517-7 




+00 06 03 6 


22 


14 ± 


07 





48 ± 


08 


4517-8 


12:32:43.44 


+00:07:08.4 


22 


20 ± 


04 


1 


22 ± 


06 


4517-9 


12:32:42.72 


+00:05:38.4 


22 


21 ±0 


02 





67 ±0 


02 


4517-10 


12:32:46.32 


+00:06:28.8 


22 


24 ±0 


03 





86it0 


03 


4517-11 


12:32:41.52 


+00:06:25.2 


22 


33 ± 


21 





52 it 


22 


4517-12 


12:32:44.16 


+00-06-32 4 


22 


43 ± 


11 


1 


42 it 


13 


4517-13 


12:32:43.92 


+00-06-32 4 


22 


56 ± 


04 


1 


26 ± 


05 


4517-14 


12:32:43.20 


+00:06:28.8 


22 


57 ±0 


09 





33 ± 


12 


4517-15 


12:32:42.96 


+00:06:21.6 


22 


89 ±0 


06 


1 


11 ±0 


06 


4517-16 


12:32:43.68 


+00:06:21.6 


23 


07 ± 


03 





98 it 


05 


4517-17 


12:32:46.08 


+00-06-28 8 


23 


08 ± 


13 





90 it 


16 


4517-18 


12:32:44.40 


+00:06:25.2 


23 


18 ± 


05 





69 it 


09 


4517-19 


12:32:41.04 


+00:04:58.8 


23 


31 ±0 


06 


1 


03 ± 


09 


4517-20 


12:32:48.48 


+00:05:27.6 


23 


38 ±0 


07 





55itO 


08 


4517-21 


M-VtM 76 


+00-05-24 


23 


50 ± 


06 


1 


93 it 


11 


4517-22 


12:32:45.12 


+00-06-25 2 


23 


62 ± 


10 





63 it 


13 


4517-23 


12:32:46.32 


+00-06-36 


23 


64 ± 


10 





34 it 


13 


4517-24 


12:32:45.36 


+00:06:36.0 


23 


71it0 


08 





56 ± 


12 


4517-25 


12:32:44.88 


+00:06:25.2 


23 


72 ±0 


07 





42 ± 


10 


4517-26 


12:32:43.20 


+00-07-08 4 


23 


73 ± 


09 





36 it 


11 


4517-27 


12:32:45.84 


+00-06-21 6 


23 


78 ± 


12 





74 it 


15 


4517-28 


12-32-45 36 


+00-06-32 4 


23 


80 ± 


09 





64 it 


11 


4517-29 


12:32:43.20 


+00:07:12.0 


23 


84 ±0 


07 





35it0 


09 


4517-30 


12:32:45.60 


+00:04:55.2 


23 


84 ±0 


25 


1 


58it0 


29 


4517-31 


12:32:45.84 


+00:06:28.8 


23 


85 ± 


10 





62 ± 


14 


4517-32 


12:32:45.36 


+00:06:32.4 


23 


86 ± 


08 





32 it 


11 


4517-33 


12-32-45 60 


+00-06-25 2 


23 


87 ± 


08 


1 


04 it 


11 


4517-34 


12:32:43.44 


+00:06:21.6 


23 


89 ±0 


12 


1 


88it0 


14 


4517-35 


12:32:43.20 


+00:06:32.4 


23 


92 ±0 


07 


1 


69it0 


08 


4517-36 


12:32:45.60 


+00:06:25.2 


23 


95 ± 


13 





36 ± 


17 


45 17-37 


12:32:44.40 


+00-06-25 2 


24 


02 ± 


08 





38 it 


14 


4517-38 


12:32:42.00 


+00-06-25 2 


24 


07 ± 


16 


1 


48 it 


18 


4*^17 30 


19*39'4'^ 84 


-l-^ft•^^^•'^^^ ft 


24 


ns -t- n 
uo in u 


1^1 


n 
u 


70 -1- n 
1 u m u 




4517^0 


12:32:44.40 


+00:06:21.6 


24 


10 ±0 


10 





58it0 


13 


4517^1 


12:32:42.96 


+00:06:14.4 


24 


19 ±0 


19 





82 ± 


21 


4517^2 


12:32:44.40 


+00:06:10.8 


24 


20 ±0 


11 





90it0 


15 


4517^3 


12:32:43.44 


+00:06:25.2 


24 


23 ±0 


19 





65it0 


22 


4517^ 


12:32:45.84 


+00:06:36.0 


24 


24 ±0 


11 





95it0 


14 


4517^5 


12:32:45.84 


+00:06:32.4 


24 


28 ±0 


08 





73it0 


13 


4517^6 


12:32:45.84 


+00:06:36.0 


24 


30 ±0 


16 





49 ± 


17 


4517-47 


12:32:43.20 


+00:06:32.4 


24 


31 ±0 


11 





52 ± 


16 


4517^8 


12:32:44.40 


+00:06:25.2 


24 


35 ±0 


19 





53 ± 


21 


4517-49 


12:32:45.60 


+00:06:32.4 


24 


41 ±0 


12 





46it0 


18 


4517-50 


12:32:45.36 


+00:06:28.8 


24 


45 ±0 


11 





80it0 


13 
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Table B4. Photometry and astrometry of the 50 brightest globular cluster 
candidates on the WFPC2 frames of IC 5 176. The object Ust is sorted on V 
magnitude (brightest first). 



ID 


RA 
(J2000) 


DEC 
(J2000) 




V 
mag 






V-I 
mag 




5176-1 


22 


14:54.96 


-66:50:38.4 


21 


97 ±0 


04 


1 


05 ±0 


04 


5176-2 


22 


14:43.44 


-66:51:03.6 


22 


18 ±0 


06 





90 ±0 


07 


5176-3 


22 


1 s-no no 


-66:51:14.4 


22 


41 ±0 


02 


1 


82 ±0 


02 


5176-4 


22 


15:05.04 


-66:50:09.6 


2.3 


25 ±0 


04 


1 


18 ±0 


05 


5176-5 


22 


14:51.84 


-66:51:18.0 


23 


25 ±0 


04 


1 


13 ±0 


04 


5176-6 


22 


15:01.68 


-66:50:34.8 


23 


33 ±0 


06 


1 


19 ±0 


08 


5176-7 


22 


14:56.88 


-66:51:25.2 


23 


54 ±0 


05 


1 


17 ±0 


05 


5176-8 


22 


1 5 on 96 


-66:50:42.0 


23 


63 ±0 


04 


1 


04 ±0 


07 


5176-9 


22 


14:52.08 


-66:51:00.0 


23 


68 ±0 


04 


1 


02 ±0 


07 


5176-10 


22 


14:51.12 


-66:51:324 


23 


76 ±0 


06 


1 


14 ±0 


07 


5176-11 


22 


14:50.16 


-66:51:21.6 


23 


79 ±0 


04 


1 


07 ±0 


06 


5176-12 


22 


14:51.36 


-66:51:21.6 


23 


81 ±0 


05 


1 


01 ±0 


06 


5176-13 


22 


14:57.12 


-66:51:18.0 


23 


81 ±0 


04 


1 


17 ±0 


05 


5176-14 


22 


14:54.96 


-66:51:324 


23 


82 ±0 


05 


1 


42 ±0 


06 


5176-15 


22 


15:02.88 


-66:49:51.6 


23 


91 ±0 


08 





73 ±0 


10 


5176-16 


22 


1 A -CO r\o 

14:58.08 


-66:51:39.6 


23 


93 ±0 


11 





92 ±0 


11 


5176-17 


22 


14:59.76 


-66:49:55.2 


24 


00 ±0 


09 





84 ±0 


10 


5176-18 


22 


15:02.64 


-66:50:09.6 


24 


20 ±0 


06 





98 ±0 


07 


5176-19 


22 


14:52.80 


-66:51:00.0 


24 


27 ±0 


06 


1 


11 ±0 


07 


5176-20 


22 


14:53.52 


-66:50:56.4 


24 


48 ±0 


11 


1 


30 ±0 


12 


5176-21 


22 


14:58.56 


-66:49:55.2 


24 


49 ±0 


06 





83 ±0 


08 


5176-22 


22 


14:56.64 


-66:51:18.0 


24 


53 ±0 


09 


1 


09 ±0 


10 


5176-23 


22 


14:51.36 


-66:51:18.0 


24 


57 ±0 


07 





89 ±0 


10 


5176-24 


22 


14:51.60 


-66:51:07.2 


24 


60 ±0 


09 


1 


05 ±0 


11 


5176-25 


22 


14:56.64 


-66:51:39.6 


24 


66 ±0 


08 





78 ±0 


09 


5176-26 


22 


15:03.84 


-66:50:20.4 


24 


67 ±0 


06 


1 


11 ±0 


08 


5176-27 


22 


14:53.52 


-66:50:56.4 


24 


71 ±0 


10 





40 ±0 


16 


5176-28 


22 


14:51.84 


-66:51:18.0 


24 


71 ±0 


09 





58 ±0 


11 


5176-29 


22 


15:00.24 


-66:50:384 


24 


74 ±0 


10 


1 


14 ±0 


11 


5176-30 


22 


14:56.40 


-66:51:25.2 


24 


81 ±0 


12 





90 ±0 


13 


5176-31 


22 


14:54.72 


-66:50:42.0 


24 


83 ±0 


11 





87 ±0 


13 


5176-32 


22 


14:56.64 


-66:50:24.0 


24 


86 ±0 


12 





95 ±0 


14 


5176-33 


22 


15:01.92 


-66:50:13.2 


24 


92 ±0 


07 


1 


28 ±0 


11 


5176-34 


22 


14:55.92 


-66:50:24.0 


24 


92 ±0 


09 





90 ±0 


11 


5176-35 


22 


14:53.04 


-66:50:52.8 


24 


93 ±0 


14 





46 ±0 


24 


5176-36 


22 


14:51.12 


-66:51:18.0 


24 


97 ±0 


08 





75 ± 


09 


5176-37 


22 


15:05.04 


-66:50:09.6 


25 


03 ±0 


10 


1 


40 ±0 


11 


5176-38 


22 


14:55.68 


-66:50:31.2 


25 


09 ±0 


20 


1 


30 ±0 


22 


5176-39 


22 


14:55.68 


-66:50:31.2 


25 


09 ±0 


20 


1 


29 ±0 


22 


5176-40 


22 


14:51.84 


-66:51:21.6 


25 


09 ±0 


11 


1 


15 ±0 


12 


5176^1 


22 


15:04.08 


-66:49:48.0 


25 


10 ±0 


09 





37 ±0 


12 


5176-42 


22 


14:54.96 


-66:51:36.0 


25 


18 ±0 


09 


1 


37 ±0 


11 


5176-43 


22 


14:50.64 


-66:51:18.0 


25 


21 ±0 


11 





78 ±0 


13 


5176-44 


22 


15:01.68 


-66:50:45.6 


25 


30 ±0 


09 


1 


13 ±0 


11 


5176-45 


22 


15:01.20 


-66:50:45.6 


25 


32 ±0 


15 


1 


09 ±0 


18 


5176^6 


22 


14:59.52 


-66:49:58.8 


25 


35 ±0 


10 





96 ±0 


14 


5176-47 


22 


15:01.68 


-66:50:56.4 


25 


42 ±0 


17 


1 


28 ±0 


19 


5176-48 


22 


14:52.56 


-66:50:13.2 


25 


42 ±0 


19 


1 


32 ±0 


19 


5176-49 


22 


14:51.36 


-66:51:21.6 


25 


46 ±0 


11 


1 


51 ±0 


13 


5176-50 


22 


14:54.72 


-66:50:20.4 


25 


49 ±0 


09 


1 


09 ±0 


10 



Table B5. Photometry and astrometry of the globular cluster candidates on 
the WFPC2 frames of NGC 7814. The object list is sorted on V magnitude 
(brightest first). 



ID 


RA 
(J2000) 


DEC 
(J2000) 




V 
mag 








V-I 
mag 




7814-1 


00:03 


13.20 


+16:08:27.6 


21 


19 ± 





03 





78 ±0 


08 


7814-2 


00:03 


09.84 


+16:06:50.4 


21 


34 ± 





03 





61 ±0 


03 


7814-3 


00:03 


10.80 


+16:08:02.4 


21 


43 ± 





15 





43 ±0 


17 


7814-4 


00:03 


08.64 


+ 16:08:49.2 


21 


54 ± 





10 





68 ±0 


10 


7814-5 


00:03 


13.44 


+16:08:24.0 


22 


15 ± 





08 


1 


12 ±0 


09 


7814-6 


00:03 


12.72 


+16:08:204 


22 


44 ± 





07 





60 ±0 


12 


7814-7 


00:03 


11.28 


+16:07:55.2 


22 


67 ± 





17 





80 ±0 


18 


7814-8 


00:03 


12.24 


+16:08:42.0 


22 


73 ± 





06 





30 ±0 


13 


7814-9 


00:03 


12.48 


+16:07:26.4 


23 


03 ± 





04 


1 


05 ±0 


06 


7814-10 


00:03 


10.80 


+16:08:45.6 


23 


09 ± 





19 





84 ±0 


22 


7814-11 


00:03 


14.16 


+16:06:504 


23 


27 ± 





10 





88 ±0 


11 


7814-12 


00:03 


13.20 


+16:06:57.6 


23 


27 ± 





03 


1 


09 ±0 


04 


7814-13 


00:03 


15.12 


+16:08:13.2 


23 


32 ± 





08 


1 


05 ±0 


09 


7814-14 


00:03 


12.96 


+16:07:30.0 


23 


42 ± 





03 


1 


72 ±0 


14 


7814-15 


00:03 


13.68 


+16:08:24.0 


24 


12 ± 





14 


1 


17 ±0 


20 


7814-16 


00:03 


13.44 


+16:08:24.0 


24 


32 ± 





11 





76 ±0 


16 


7814-17 


00:03 


13.44 


+16:08:34.8 


24 


63 ± 





20 


1 


92 ±0 


22 
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